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ABSTRACT We study two different methods to test Gaussianity in CMB maps. One of them 
is based on the partition function and the other on the morphology of hot and cold spots. The 
partition function contains information on all the moments and scales, being a useful quantity to 
compress the large data sets expected from future space missions like Planck. In particular, it 
contains much richer information than the one available through the radiation power spectrum. 
The second method utilizes morphological properties of hot and cold spots such as the eccentricity 
and number of spots in CMB maps. We study the performance of both methods in detecting 
non-Gaussian features in small scale CMB simulated maps as those which will be provided by the 
Planck mission. 
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1. INTRODUCTION 

Future CMB experiments like the Planck mission will provide very large data sets at 
small angular scales with a high sensitivity. Their analysis will require the develop- 
ment of new and sophisticated methods capable of managing the data, performing 
an optimal separation of the different components (CMB, Galactic and extragalac- 
tic foregrounds), extracting their statistical properties and finally determining the 
fundamental cosmological parameters. Two methods for foreground removal have 
already been proposed based on Maximum Entropy (Hobson et al. 1998, 1999a) 
and Wiener filter (Tegmark and Efstathiou 1996). The power spectrum of the tem- 
perature fluctuations is the most widely used quantity in order to compress the data 
and obtain the cosmological parameters. This quantity only carries information on 
the second moment of the distribution and completely characterizes the statistical 
properties of the data if these are Gaussian distributed. However, at present there 
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is no experimental evidence for that statistical behaviour in the small scale tem- 
perature fluctuations. Moreover, topological defect models predict non-Gaussian 
temperature fluctuations, ft is therefore very important to explore other quanti- 
ties that can be sensitive to non-Gaussian features and can test the assumption of 
Gaussianity usually made for the cosmological signal. 

Several studies of the Gaussian character of the CMB signal based on different 
statistical and morphological quantities have already been performed. They include 
extrema correlation function (Kogut et al. 1996, Barreiro et al. 1998), number of 
hot and cold spots (Coles and Barrow 1987), genus (Gott et at. 1990, Kogut et 
al. 1996), bispectrum (Ferreira et al. 1998, Heavens 1998), wavelets (Hobson et al. 
1999b, Pando et al. 1999) and multifractals (Pompilio et al. 1995). 

In this paper we concentrate on the recently proposed partition function of the 
CMB temperature and also on the morphology of the temperature extrema. The 
usefulness of the partition function to study the Gaussianity of the CMB has already 
been pointed out in a previous work (Diego et al. 1999). The study of morphological 
properties of CMB extrema, such as number, curvature and eccentricity, has been 
done by Barreiro et al. (1997) for Gaussian temperature fluctuations. Here we 
will probe the practical performance of both methods by testing their abilility to 
find the non-Gaussian features of very weakly perturbed Gaussian random fields. 
In particular we will simulate small scale temperature maps by slightly perturbing 
Gaussian ones and generating a small: a) skewness or b) kurtosis. 

2. TESTS OF GAUSSIANITY 

2.1 The partition function 

The partition function has recently been introduced to study CMB data by 
Diego et al. (1999). Below we concentrate on the application of this function to the 
study of Gaussianity. 

The partition function is defined as follows: 



The quantity Hi(5) is called the measure, it is a function of 5 which is the size 
or scale of the boxes used to cover the sample. The boxes are labeled by i and 
^boxcs(<5) is the number of boxes (or cells) needed to cover the map when the grid 
with resolution 8 is used. The exponent q is a continuous real parameter that plays 
the role of the order of the moment of the measure. 

Let us consider a CMB map of N pixels. Now the map is divided in boxes of 
size S x S pixels and the measure Hi{5) is computed in each one of the resulting 
boxes. Changing both, q and <5, one calculates the function Z(q, 5). We would like 
to emphasize that the calculation of Z(q, S) is O(N). 

One is free to make any choice of the measure /i(5) provided that several conditions 
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are satisfied, the most restrictive being ^i(S) > 0. There arc no general rules to 
decide which is the best choice. To test Gaussianity we use the following measure: 

pix^GboXi 

Thus the measure in the box i is the sum of the absolute temperatures T p ; x of the 
pixels inside the box in units of Kelvin. The constant T» is a normalization constant. 
The measures are interpreted as probabilities and they have to be normalized, i.e 
Hi = l. So T* is simply the sum of the absolute temperatures over all pixels and 
therefore is a constant for all boxes and scales. 



Summarizing, Z(q, 5) contains information at different scales and moments. The 
multi-scale information gives an idea of the correlations in the map, meanwhile the 
moments are sensitive to possible asymmetries in the data, as well as some deviations 
from Gaussianity. In order to study the last property let us introduce the generating 
function: 

Nboxcs(<5) 

G„(g) = <e«"> = m ]T e«*< 4 >. (3) 

* boxes / -j 

Defining a new function H(q) in the following way, 
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H(q) = InG^q) - q{») - S-Za, (4) 

this function vanishes for any value of q when a Gaussian distributed variable is 
considered. In the case of the CMB temperature, if the universe went through a 
phase of inflation then the fluctuations generated in the CMB will be Gaussian. 
However, due to the cosmic and sample variance H (q) for a given realization will 
not exactly vanishes, as we will see in the next section. 

2.2 The morphology of hot and cold spots 

The properties of the CMB temperature peaks for the Gaussian case have been 
recently studied in Barreiro et al. (1997). The local description of maxima in- 
volves the second derivatives of the field along the two principal directions. As 
usual, the curvature radii are defined by R x = [-A' 1 '(maa;)/2]^ 1 / 2 and R 2 = 
[— A 2 (max)/2]~ 1 / 2 , where A is the temperature field normalized to the rms fluc- 
tuations. Then, we can associate two invariant quantities to any maximum: the 
Gaussian curvature k and the eccentricity e given by 
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In addition to the Gaussian curvature and the eccentricity another interesting 
property is the number of maxima (minima) N >v above (below) a given temper- 
ature threshold v (—v) in CMB maps. For symmetric distributions, such as the 
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FIGURE 1. Top: The function H(q) for two different sizes of the boxes. Bottom: Distribution 
of Ln(Likelihood) corresponding to H(q) for the Gaussian and non-Gaussian models (see 
text). 



Gaussian one, the number of maxima and minima are statistically undistinguish- 
able. Thus, any asymmetry present would be a mark of non-Gaussianity. Instead 
of using those morphological quantities for extrema in practice it is simpler to use 
excursion sets, i.e. hot and cold spots regions above and below a given temperature 
threshold. Therefore, we will count the number of excursion sets from the maps. To 
estimate the eccentricity of the excursion sets we fit each region (formed by at least 
6 pixels) with an ellipse of equal area as the corresponding excursion set. Finally, 
the Gaussian curvature is a more difficult quantity to estimate from the simulated 
maps and we will leave it for a forthcoming paper. 



3. SIMULATIONS 

We perform small scale simulations of size 12.8 x 12.8 degrees with pixel and 
antenna sizes similar to those of the best Planck channels, i.e. 1'.5 and 5' respec- 
tively. The amount of instrumental noise in the maps provided by this mission is 
expected to be much smaller than the CMB signal, then, as a first approximation, it 
is not considered in this work. We simulate maps with Gaussian and non-Gaussian 
statistical properties, that will be used to test the previously discussed methods. 
The properties of the Gaussian perturbations are completely defined by the stan- 
dard inflationary model. Although many models have been proposed for which 
non-Gaussian temperature perturbations are produced, however all their statistical 
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properties are either not given or are very complicated to simulate. We then consider 
"generic" non-Gaussian perturbations derived from Gaussian ones using transfor- 
mations given by Cole and Weinberg (1992). The non-Gaussian temperature fields 
are derived from Gaussian fields in two ways. First, we perform the following trans- 
formation to an uncorrelated 2D Gaussian field G: exp(aG) — exp(a 2 /2) with a a 
parameter which recovers the Gaussian case when a — > 0. An uncorrelated non- 
Gaussian field with zero mean is thus obtained. The radiation power spectrum 
corresponding to a standard CDM model normalized to COBE/DMR is introduced 
by rescaling the amplitude of each Fourier mode to that spectrum. By rescaling 
the power spectrum the 1-point distribution is modified with respect to the uncor- 
related one, but it is still non-Gaussian. Second, we consider the transformation: 
G[exp(aG) + exp(-aG)]. Then the same radiation power spectrum is introduced 
as before. Finally the maps are smoothed with a FWHM= 5'. In the first case 
we construct a non-Gaussian field whose distribution is skewed with a positive tail 
(P) and in the second one the distribution is broadened with both positive and 
negative tails (B). For both non-Gaussian fields we consider small enough values 
of the parameter a such that their 1-point distributions are undistinguishable from 
a Gaussian. In addition their power spectra are undistinguishable from the one of 
the Gaussian CDM model by construction. 

4. RESULTS 

The results of comparing two sets of simulated maps, one corresponding to 
Gaussian fluctuations and the other to one of the previous non-Gaussian fields, 
using the methods described in section 2 are summarized in figures 1 and 2 for 
the functions H(q) and N >u , respectively. We have considered a value a = 0.8 for 
which it is impossible to distinguish between the corresponding 1-pdf. For the set of 
values of H (q) (N >u ) given by the top of figure 1 (2) we have constructed a likelihood 
function whose values for each of the simulated Gaussian and non-Gaussian maps 
(compared to the Gaussian model) are shown in the bottom of figure 1 (2). We see 
that a positive result, i.e. it is posible to distinguish between the two distributions 
for most of the simulations, is obtained for the P non-gaussian field for both the 
H(q) and N >v quantities. In the case of the B non-gaussian field the difference is 
less pronounced but it is still possible to discriminate between both models in a 
number of cases. 

We have also performed a similar test for the eccentricity e>„ but in this case it is 
not possible to discriminate the Gaussian and any of the non-Gaussian fields. 

Summarizing, the function H(q) as well as iV>„ are sensitive quantities with 
which to probe the possible non-Gaussian features that can be present in future 
high sensitivity maps. On the contrary, the eccentricity is not a good discriminator 
of Gaussianity at least for the small maps considered in the present work. 
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FIGURE 2. Top: N >u for thresholds i/ in the range 1 — 4. Bottom: Distribution of 

Ln(yLikelihood) corresponding to -/V>„ for the Gaussian and non-Gaussian models (see 
text). 
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